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Using extremely deep (rms ^^S.SjtiJy/bm) 1.4 GHz sub-arcsecond resolution MERLIN+VLA ra- 
dio observations of a 8^5x8^5 field centred upon the Hubble Deep Field North, in conjunction 
with Spitzer 24 jim data we present an investigation of the radio-MIR correlation at very low flux 
densities. By stacking individual sources within these data we are able to extend the MIR-radio 
correlation to the extremely faint (^microJy and even sub-microJy) radio source population. Ten- 
tatively we demonstrate a small deviation from the correlation for the faintest MIR sources. We 
suggest that this small observed change in the gradient of the correlation is the result of a suppres- 
sion of the MIR emission in faint star-forming galaxies. This deviation potentially has significant 
implications for using either the MIR or non-thermal radio emission as a star-formation tracer at 
low luminosities. 
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1. Introduction 

Radio and infrared emission from galaxies in both the nearby and distant Universe is thought 
to arise from processes related to star-formation, hence resulting in the correlation between these 
two observing bands. The infrared emission is produced from dust heated by photons from young 
stars and the radio emission predominately arises from synchrotron radiation produced by the ac- 
celeration of charged particles from supernovae explosions. It has however recently been suggested 
that at low flux density and luminosities there may be some deviation from the tight well-known 
radio-IR correlation seen for brighter galaxies [|l|, §]. 

|[T]] argue that while the IR emission from luminous galaxies will trace the majority of the 
star-formation in these sources, in low luminosity galaxies the IR emission will be less luminous 
than expected considering the rate of star-formation within the source (i.e. the IR emission will not 
fully trace the star-formation). In this scenario the reduced efficiency of IR production relative to 
the source star-formation rate (SFR) would be the result of inherently lower dust opacities in lower 
luminosity sources and consequently less efficient reprocessing of UV photons from hot young 
stars into IR emission. The simple consequence of this is that at lower luminosities the near linear 
radio-IR correlation Lradio 

ocL^j^, with y > 1 (e.g. ^) will be deviated from. Of course such an 
assertion is dependent upon the radio emission providing a reliable tracer of star-formation at low 
luminosities which may be equally invalid. 

Recently have presented a statistical analysis of Australia Telescope Compact Array (ATCA) 
20 cm observations of the 24 /xm sources within the Spitzer Wide Field Survey (SWIRE). In this 
work ^ have co-added sensitive (rms'^SO /iJy) radio data at the locations of several thousand 
24 /im sources. Using this method they have statistically detected the microJy radio counterparts 
of faint 24 /xm sources. At low flux densities (S24^m = 100 /xJy) they confirm the IR-radio cor- 
relation but find it to have a lower coefficient (Si.4ghz = 0.039 S24^m) than had previously been 
reported at higher flux densities. This coefficient is significantly different from results previously 
derived from detections of individual objects (e.g. ^) and is speculated by [Q] to be the result of a 
change in the slope of the radio-IR correlation at low flux densities. 

In this work (which is descirbed in more detail in [|6|, |7|]) we have utilised very deep, high 
resolution 20 cm observations of the Hubble Deep Field North and surrounding area made using 
MERLIN and the VLA [§] in combination with publicly available 24 /xm Spitzer source catalogues 
from GOODS to study the MIR-Radio correlation for microjansky radio sources. This study ex- 
tends the flux density limits of the radio-IR correlation by more than an order of magnitude for 
individual sources and overlaps the flux density regime studied using statistical stacking methods 
by other authors. 

2. Observations 
2.1 Radio Data 

Extremely deep radio observations of the HDF-N region were made in 1996-97 at 1.4 GHz 
using both MERLIN and the VLA. These observations were initially presented in [||], [Q] and [p^]. 
The results from the combined 1 8 day MERLIN and 42 hr VLA observations are described in detail 
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Figure 1: Mean (left) and median (right) images of the 1.4 GHz radio emission for all sources within the six 
faintest 24 jim flux density logarithmic bins plotted in Fig. 1 in descending flux density order from top-left to 
bottom-right. The range of 24 /xm flux density over which each image has been stacked is shown at the top 
of individual panels. Each image is contoured with levels of —2, —1.414, —1,1, 1.414, 2, 2.828, 4, 5.657, 8, 
11.31, 16, 22.63 and 32 times 3x(3.3/VN) ^Jybm"', where N equals the number of 24 /im source positions 
averaged in the map. The peak flux density, lowest plotted contour and number of IR sources which have 
been averaged over (N) is Usted at the bottom of each image panel. 



in [^. The combined MERLIN+VLA image has an rms noise level of 3.3 /iJy per O'.'l circular 
beam making it amongst the most sensitive, high-resolution radio images made to date. 

2.2 GOODS-N Spitzer 24 jUm observations 

As part of the GOODS enhanced data release^ (DRl-i- February 2005) a catalogue of Spitzer 
24 jUm source positions and flux densities for the GOODS-N field were released. This source 
catalogue is limited to flux densities >80 jxly providing a highly complete and reliable sample. At 
the time of writing this 24 jxra source catalogue represents the most complete and accurate mid- 
infrared source list publicly available for the GOODS -N/HDF-N field. 

All 24 /im sources which are detected optically in GOODS HST ACS images show an accurate 
astrometric alignment with their optical counterparts implying that the astrometry between these 
two data-sets and their subsequent catalogues is self-consistent. However, a comparison of the 
astrometric alignment of the positions of sources catalogued by GOODS derived from their HST 
ACS images [ [IT| , ^ shows there to be a systematic offset in declination of —0^342 from the radio 
reference frame. This linear declination correction, although small relative to the Spitzer resolution 

'http://www.stsci.edu/science/goods/DataProducts/ 
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at 24 /im is significant when compared to these high resolution radio data. This linear correction 
has been applied to the Spitzer source positions prior to all comparisons between the two data sets. 




Figure 2: Left: Radio 1.4 GHz versus the 24 flux density. Sources from the 8^5x8 '5 HDF-N field are 
plotted individually (small black triangles). The median radio flux density logarithmically binned by 24 /im 
flux densities are plotted as large filled red circles. The solid red line represents the best-fit line to the binned 
HDF-N data alone. Sources within the CDFS-SWIRE field detected at both 24 jim and 1.4 GHz from Norris 
et al. are plotted as either green open stars (identified as AGN) or green triangles (not identified as 
AGN). All sources detected at both 1.4 GHz and 24 /xm in the Spitzer First Look Survey (FLS) with source 
position separations of < 1 "5 are plotted in orange Jl^ ] . Note the quantisation of the SWIRE and FLS points, 
in this and subsequent plots is a result of the accuracy of the flux densities tabulated in the literature. The 
blue pluses and fitted lines in the low portion of the plot show the IR-radio correlation derived from stacking 
radio emission at the positions of 24 fim sources in the CDFS and ELAIS field by Boyle et al. ||^]. 
Right: In the upper panel the flux density ratio (S24/im/Si.4GHz) versus 24 jj.m flux density is shown. The 
individual sources from the HDF-N field are plotted as small black triangles, the median values of these 
HDF-N source binned as a function of S24/jm as filled circles, green stars show the median binned values 
for sources listed as non-AGN within the CDFS-SWIRE sample of Norris et al. [ p^ and the blue 'crosses' 
and 'pluses' show the flux density ratios derived from the stacking analysis of the SWIRE-CDFS and ELAIS 
fields respectively from Boyle et al. [^] . The overlaid black dotted line is the mean value for (S24/im/S i.4ghz) 
derived by Appleton et al. This line is equivalent to q24=0.84±0.28 with the gray filled box representing 
the area enclosed by these errors and the flux density range investigated by Appleton et al.. The values of q24 
against 24 jj.m flux density are plotted in the lower panel. The symbols within this plot are identical to those 
in the upper panel, individual HDF-N sources are not included for clarity. The additional diagonal solid 
green line represents a line of constant radio flux density of 100 /xJy, the lowest flux density of sources in 
the SWIRE sample plotted here (sources above this line are excluded by this limit from the SWIRE sample). 
This flux density limit will significantly effect the values of the binned SWIRE data points (green crosses) 
negatively biasing the four lowest flux density bins. This bias only effects the SWIRE sample. 
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3. Results & Discussion 

Using these two highly sensitive data sets the radio emission from a sample of faint 24 ;Um 
galaxies has been invesitigated. Of the 377 Spitzer sources within the radio field 303 were found to 
have total radio flux denisties in excess of 3 times the local rms in our deep radio imaging. Many 
of these sources, however, have peak flux densities comparable or below the point source detection 
threshold of these radio data. 

By using both statistically stacked images and radio flux density measurements at the positions 
of individual Spitzer 24 /im sources we have investigated the radio structures and flux densities of 
faint IR sources. In figure |l| the statistically averaged radio emission from faint 24 /im galaxies in 
the GOODS-N field are shown. Each of these images is the statistically combined radio emission 
from the location of several tens of Spitzer 24 /im sources and has been contoured at multiples of 
three times 3.3 /xJy beam^' divided by VN where N equals the number of Spitzer source positions 
that have been stacked together. As can be seen in these images the off-source noise levels achieved 
approaches the value expected when co-adding multiple images with near-Gaussian noise proper- 
ties. The co-added image rms achieved in the faintest 24 /im flux density bin (80.0 to 100.25 /iJy) 
is 0.45 and 0.56 ^uJybeam^^ in the mean and median co-added images respectively. 

The Gaussian fitted sizes of the radio emission in the stacked images (Fig|T]) provide an upper 
limit on the average size of the radio counterparts of these faint IR sources. The largest angular 
sizes of the radio emission in the median stacked images created from this sample range between 
I'.' 4 and 2". This is approximately equivalent to a linear size of lOkpc at redshifts beyond 1. These 
upper limits on the radio source sizes are consistent with radio emission on galactic and sub-galactic 
scales and originating within kpc-scale starburst systems. 

By combining the Spitzer 24 /im flux densities and the extracted 1.4 GHz radio flux densities 
for these sources it is possible to begin to characterise the radio emission of the faintest IR galax- 
ies. Shown in the left-hand portion of figure|2| is 1.4 GHz total flux density of sources within the 
8f5 X 8f5 MERLIN+VLA field plotted against their 24 /xm flux density. As a direct comparison also 
overlaid on this diagram are source flux densities from various other deep multi-wavelength obser- 
vational campaigns of different regions of sky. As can be seen these new HDF-N/GOODS-N data 
significantly extend the observed IR-radio correlation for galaxies between these two wavelengths. 

Binning and re-plotting these data in terms of S24/jm/Si.4GHz and q24 (log(S24/im/Si.4GHz) (see 
figure^ right) clearly show the slope and tightness of this correlation. However, it can be seen from 
figure^ that these data appear to show that this correlation begins to deviate for the faintest sources. 
This deviation is small (and tentative) but implies that the faintest galaxies are under-luminous 
at MIR wavelengths relative to their observed radio wavelengths. This deviate, whilst small, has 
potentially direct implications on the use of either of these bands to quantify star-formation rates 
in faint galaxies. A more complete discussion of these results and their implications can be found 
in®g]. 

4. Conclusion 

Using one of the deepest high-resolution 1.4 GHz observations made to date, in conjunction 
with deep 24 jUm Spitzer source catalogues from GOODS, we have investigated the microJy radio 
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counterparts of faint MIR sources. These observations confirm that the microJy radio source popu- 
lation follow the MIR-radio correlation and extend this correlation by several orders of magnitude 
to very low flux densities and luminosities, and out to moderate redshifts. This extension of the 
MIR-radio correlation confirms that the majority of these extremely faint radio and 24 /im sources 
are predominantly powered by star-formation with little AGN contamination. 

Statistically stacking the radio emission from many tens of faint 24 /im sources has been used 
to characterise the size and nature of the radio emission from very faint IR galaxies well below the 
nominal radio sensitivity of these data. Using these methods the MIR-radio correlation has been 
further extended and a tentative deviation in this correlation at very low 24 flux densities has 
been identified. 
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